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Abstract. We investigate the relation between the diagonal (a xx ) and off- 
diagonal (cr xy ) components of the conductivity tensor in the quantum Hall system. 
We calculate the conductivity components for a short-range impurity potential 
using the linear response theory, employing an approximation that simply replaces 
the self-energy by a constant value —ih/(2r) with t the scattering time. The 
approximation is equivalent to assuming that the broadening of a Landau level 
due to disorder is represented by a Lorentzian with the width T = H/(2t). 
Analytic formulas are obtained for both a xx and a xy within the framework of 
this simple approximation at low temperatures. By examining the leading terms 
in tr xx and a xy , we find a proportional relation between da xy /dB and Ba xx . 
The relation, after slight modification to account for the long-range nature of the 
impurity potential, is shown to be in quantitative agreement with experimental 
results obtained in the GaAs/AlGaAs two-dimensional electron system at the low 
magnetic-field regime where spin splitting is negligibly small. 
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1. Introduction 



The experimental finding by Chang and Tsui [l| of the striking similarity between 
the longitudinal resistivity p xx and the derivative of the Hall resistivity with respect 
to the electron density n e , dp xy /dn e , in the quantum Hall regime has attracted 
considerable interest and has since been a subject of a number of experimental 
0, H 0, HI IS 0, H and theoretical [S ES E3 studies. Using a low-carrier-density 



(n e <lxl0 -15 m~ 2 ) high-mobility (/i > 300 m 2 V _1 s _1 ) two-dimensional electron 
system (2DES) in GaAs/AlGaAs, Stormer et al. 0] showed that all features in p xx 
(including overshooting flanks around quantum Hall states) are faithfully reproduced 
by the derivative of p xy with respect to the magnetic field B in the form 

B^^Pp xx , (1) 

where (3 is a sample-dependent constant value (typically between 20 and 40). Note, 
as pointed out in Ref. |l|, that the differentiation by B and that by n e are basically 
equivalent to each other, — B(d/dB) = n e (d/dn e ), if the relevant variable in the 
problem is the filling factor v = n e h/eB and not n e or B separately. 

The origin of the intriguing empirical relation Eq. fl]) remains largely enigmatic. 
A possible explanation is given by Simon and Halperin [l0( , who ascribed the relation 
to the microscopic inhomogeneity in the electron density n e inevitably present in real 
2DES samples. Noting that the macroscopic value of p xx measured in experiments 
is mainly determined by the fluctuation in the local Hall resistivity p xy (r) resulting 
from the inhomogeneity in n e rather than by the local longitudinal resistivity p xx (f), 
their theory leads to Eq. (fl]) for not too low temperatures if disorders are taken into 
consideration on multiple length scales. In a recent experiment by Pan et al. Q using 
an ultrahigh mobility (p = 3100 m 2 V _1 s —1 ) 2DES at an extremely low temperature 
(~ 6 mK), experimentally measured value of p xx was interpreted [8,El| as essentially 



reflecting the difference in p xy (r) between the voltage probes placed under slightly 
different (~0.5 %) electron density n e , and accordingly as virtually irrelevant to the 
local resistivity p xx (r). Note, however, that the van der Pauw geometry used in their 
study is not necessarily an ideal setup for the measurement of the resistivity. 

In the present paper, we explicitly calculate the diagonal (a xx ) and the off- 
diagonal (cr xy ) components of the conductivity tensor in the quantum Hall system 
by employing the linear response theory. Although there already exist a number of 
sophisticated theories devoted to the calculation of a xx and a xy in the quantum Hall 
system (see, e.g., Refs. QESEIE!), they have not been app ied to the knowledge of 



the present authors, to the interpretation of the relation between the two components 
of the conductivity tensor exemplified by Eq. fl} . We take the effect of disorder into 
account by simply assuming the Lorentzian broadening of the Landau levels with 
the width T independent of B; this can readily be done by substituting a constant 
value — iT for the self-energy in the Green's function. Although this appears to be 
somewhat an oversimplified approximation, the Lorentzian with the B-independent 
width is suggested by a number of experiments to be a function that describes quite 
well the broadening of Landau levels due to disorder HSE3) 18, H)]. By contrast, the 



well-known self-consistent Born approximation [12j yields a semi-elliptical broadening, 
which is by far a less accurate representation of the experimentally observed Landau 
levels. A great advantage of the simple approximation employed in the present study 
is that it allows us to deduce analytic formulas for both a xx and a xy for low enough 
temperatures k^T <C £f with ep the Fermi energy. The analytic formulas, in turn, 
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provide us with a transparent way to examine the underlying relation between the 
two components. By picking out the most significant terms at high magnetic fields in 
the formulas, we find the relation 

^ - ABo*,, (2) 

with the coefficient A determined by scattering parameters and Ef [see Eq. ([35]) below 
for details]. The relation is analogous to Eq. ((T|) but with a notable difference that 
a xx enters the equation in squared form. Note that Eq. fl} can be rewritten as 

B^^(3a xx (3) 

by using the approximate relations for not too-small magnetic fields, p xy ~ l/a^j, 
and p xx « <J xx /a xy . In contrast to the previous study [10] . we have not introduced 
inhomogeneity in n e in our calculation. 

The relation between a xx and a xy found in the present study is compared with 
experimental results obtained in a GaAs/AlGaAs 2DES using the Hall-bar geometry, 
a geometry well-suited to the measurement of the resistivity. Care should be taken in 
the comparison, since our theoretical calculation is based on the short-range impurity 
potential, while the dominant scattering in a GaAs/AlGaAs 2DES is known to be 
of long-ranged. We find that Eq. ([53]) below obtained by modifying Eq. {2} to 
accommodate the long-range potential describes the experimental results remarkably 
well for the low magnetic field range where the spin splitting, the localization, the 
formation of edge states, and the electron-electron interaction can be neglected. 

The paper is organized as follows. In Sec. [2[ we introduce the Green's function 
to be employed in the later calculations. Components of the conductivity tensor are 
calculated in Sec. which are shown in |Appendix A| to approach the semiclassical 
formulas asymptotically for B — > 0. The relation between a xx and a xy is examined 
in Sec. UJ and is compared with experimental results in Sec. [5J after modification 
to account for the long-range nature of the impurity potential. The validity of our 
approximation and the magnetic-field range for our approximation to be accurate are 
discussed in Sec. El followed by concluding remarks in Sec. [7] 

2. Impurity scattering in the quantum Hall system 

We consider a 2DES in a magnetic field perpendicular to the 2D plane. The 
Hamiltonian of the system is given by 

#qh = H + V imp , (4) 

#o =-^(p + eA) 2 1 (5) 
Zm* 

where p denotes the momentum operator, — e is the charge of an electron, A is 
the vector potential of the magnetic field (0, 0, B) and Vimp represents the impurity 
potential. We neglect spins for simplicity. The term Ho in the Hamiltonian gives the 
Landau levels. The eigenfunction of Ho in the Landau gauge is given by 

<t>kN{x, y) = -^=e lkx XN(y - Vk), (6) 

where L is the length of the system, xn denotes the eigenfunction of the harmonic 
oscillator in the Nth Landau level whose energy is given by En = hu> c (N + 1/2) with 
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uj c = e\B\/m* the cyclotron frequency, and y k — —k£ 2 is the guiding center with 
I = yj h/e\B\ the magnetic length. 

We consider a short-range potential of the form 

vw(0 = $3w-3)- ( 7 ) 



Owing to the impurity potential, Landau levels acquire width, which are otherwise 
delta functions placed at e = Em (N = 0, 1, 2,...). The resulting density of states 
(DOS), or the line shape of the impurity-broadened Landau levels, has been calculated 
for various types of impurity potential. For a white-noise potential (impurities with 
constant strength Vi distributed at random positions fp, the broadening was shown 



to be well described by a Gaussian line shape [20 . l2l[ |22j . Calculations were also 
done assuming a distribution P{Vi) in the strength of the impurity scattering Vi 
HI, HI [il H, 0, 0, Hz|. Brezin et al. [HI and Benedict et al. showed 



that a Lorentzian distribution of P{Vi) results in DOS described by a Lorentzian line 
shape. Lorentzian broadening of the Landau levels is consistent with experiments on 
the tunneling into a 2DES [lg, [13 or measurement of the magnetization in a 2DES 

13, Ell- 
in the present paper, we start by assuming the Lorentzian DOS 

J w-5g£; ( .-^'). + i. - < 8 > 

As will be shown, this simple approximation allows us to deduce analytic formulas 
of the conductivity tensor, which proves to be essential for the later analysis of the 
relation between the components of the conductivity tensor. 

The simple DOS Eq. JSj) implies analogous simplicity in the electron Green's 
function. For sufficiently short-ranged impurity potential, the Green's function can 
be written in the diagonalized form as 

1 



G N (s)S N>N >5 k>k > = (N,k 



QH 



N',k')= Sn - n ' 6 ^' ( 9 ) 



where \N,k) represents the eigenstate of the unperturbed Hamiltonian given by 
Eq. ([6]), and Ejv(e) denotes the self-energy resulting from VJmp- The DOS is related 
to the imaginary part of the electron Green's function ^ by 

^) = -iE^^ ( 10 ) 

with pjv(e) introduced as 

PN (e) = -ImG N (e + iO). (11) 

It is easy to see that Eq. (fTU)) reproduces Eq. © if the self-energy T, N (e) in Eq. © is 
replaced by a constant value — iT = —ih/(2r), yielding 

G N (e + t 0) = J— -. (12) 

e — En + il 

We exploit the simple Green's function Eq. (fT2| in the following calculations. 
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Figure 1. The diagram for the current-current correlation function. 



3. Conductivity tensor 



We introduce the particle-current operator j of the form 
3 = —:(p + el). 



(13) 



The conductivity tensor a a f3 (with a and (3 representing either of x or y) of the 2DES 
is given by the Kubo formula 



o- a ff(uj) = Re 



-(K af3 (u + iO)-K a p(0)) 



(14) 



where K a p represents the thermal Green's function corresponding to the current- 
current correlation function 

„2 r h/k B T 

Are^ T {T Tja (T) ]f3 {Q)), (15) 



e 

mi 



with L the system size, T T the chronological operator and uj n = 2mrkBT/h for an 
integer n. The bracket (...) here denotes the ensemble average. In the calculation of 
the conductivity tensor (|14p . we consider only the loop diagram shown in Fig. [T] and 
neglect the correction from the current vertex part. The correlation function K a p is 
then written as 

K af3 (iuj n ) = -^j^Y. E (N,k\ Ja \N',k')(N',k'\ Jp \N,k) 

ui m N,k,N'.k' 



xGN>(ifv^m + ifuvn + £F)Gjv(ifiu; m + £p), 



(16) 

where the electron Green's function Gat is given by Eq. (fT2|) and the matrix elements 
of the particle current are 



(N,k\ Jx \N>,k>) = (-±J^±1 Sni<n+1 



(N,k\jy\N',k r 



mi 



N + l. . h N \ 

-^-<W+i + i—d-5 N , tN ^ 5 k , k , 



(17) 



Performing analytic continuation of iu n to ui and taking the limit lo — > + i 0, we 
obtain the dc parts of the diagonal and off-diagonal components in the conductivity 
tensor in the forms 



<r xx (T,s F ) = |!(^ c ) 2 |°° d £ (-^) jT(N+l)p N (s)p N+1 (e) 



(18) 



cr xy (T,e F ) = ~2nh^ c 



,. , 8G N+1 (e + iO) dG N (e + iO) 

def{e)(N + 1) p N (e) — Pn+i{s) 



(19) 

where /(e) = l/{exp[(e— e F ) / (k-BT)] + l} is the Fermi distribution function. Equations 
(I18|) and (fT9|) are basically equivalent to Eqs. (50) and (51) in Ref. [28j by Jonson and 
Girvin, except that the self-energy in the Green's function is replaced by a constant 
value in our case. 

We can calculate the components of the conductivity tensor Eqs. (|18p and (fT9|) 
further using the electron Green's function (fT2f . We first examine the diagonal 
component a xx . For k^T <C £p, we can approximate the derivative of the Fermi 
distribution function by the delta function —df(e)/de ~ S(e — e F ). Thus Eq. (fT8| 
becomes 

2 00 

^x(ef) - |^(^) 2 + 1W(£fW+i(£f). (20) 

Introducing dimensionless parameters 

F huj r 2 ' ^ ^ 



and 



7= A (22) 



we can rewrite Eq. (|2"D|) as 

e 2 7 2 ^ 2X F + 1 
^epj _ ^ i + ^ 2 (Xp _ + ^2 ■ (^) 

To evaluate the summation over N in Eq. (|23|) . we use the Poisson sum formula 
V = ±2™ Y gf TV \ e -a*»x*-**\uh 

N=-oo r ~ 1 v=-oa 



7r[— sin(27rXF) ± i sinh(27T7)] 
cosh(27T7) — cos(27rXp) ' 
where #(£) is the unit step function. Using Eq. ([24]), we derive from Eq. ([23] 



(24) 



( e \ e 2 2 7 (X F + 1/2) sinh(27r 7 ) e 2 
ft 1 + 4 7 ^ cosh(27T7j — cos(27rAp) h 

where we approximated X)w=o by StvL-oo' n °thig that terms with N < are 
negligibly small at £f in the typical situations e F T. Equation (j25|) bears the same 
form as Eq. (2.11) in Ref. [12J by Ando, if we replace our X F and 7 with X' /(huj c ) 
and X" I (Huj c ) , respectively. In the second equality in Eq. (|2"5|) , we introduced the 
notation d a p for the conductivity a a p normalized by e 2 /h. 

Next we examine the off-diagonal component a xy of the conductivity tensor. 
Introducing a variable of integration 

x = h-\ < 26 > 



and performing the integration by parts, we rewrite Eq. (119|) as 

1 x / v ' df(hw c (x + i/2)) 



p 2 j=U r°° ( , 

^S V ) = -—Y J {N + 1) JX(-- 



dX 



L N (X), 



(27) 



x r 



with 



Ljv(X) = (hw c f 



For k^T <C e F , we obtain 



flUJr X 



PN 



pn+1 hu c x 



, 1 \ \ dG N+1 (huj c (X / + 1/2) + iO) 



dX' 



, t 1\\ dG N (hm c (X' + 1/2) + zO) 
dX' 



AX'. 



^2 00 



5^ (JV + I)Ljv(Xf). 



(28) 



(29) 



N=0 



Using Eq. (p~2|) and performing the integration L(X) in Eq. (|28|) up to X F , we obtain 
2 7 3 ^ [ 1 (X F -N) 2N + 1 



Ca;y(eF) — — 



27r 2 ft 1 + 4 7 : 



E 

N=0 



E 

iV=0 



arctan 



2 7 2 (X F - N) 2 + 7 2 {X F - TV) 2 + 7 2 

X F -N\ n S 
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(30) 



We can evaluate the summation over N in Eq. (|30[) . following a similar procedure as 
in the calculation from Eqs. (f2"3")) to ([23)) . The first line in the right hand side (r.h.s.) 
of Eq. (13TJ)) becomes 



h cosh(27T7) — cos(27rXp) 



-7 sin(27rXf 



4 7 2 (X F + l/2) . 



sinh(27i"7) 



(31) 



1 + 4 7 2 

where we used Eq. (|24|) . employing the approximation X)jv=o — * Siv=-oo as before 
Along the same line, we can accurately approximate the last term in the r.h.s. of Eq 
(HB by 

9 -1 r °° 

1 x - 

> arctan 



h it 



.N=-c 



X V -N\ 7T 



(32) 



noting that arctan((X F - AO/7) - tt/2 for N < since (X F - N)/i = [sy 
(N + 1/2)Hlj c ]/T > for e F > T. Using the relation arctan((X - N)/f) 
J j/[(X - N) 2 + 7 2 ]dX + const, and Eq. (|2l|). we can rewrite Eq. fl33J) further as 

00 f x F 1 



E 



7 



o 7T (X - N) 2 + 7 2 



rdX 



arctan (coth(7T7) tan(7rX F )) + Int X F 



with Int(£) representing the integer part of £. We finally arrive at 



1 \ 1 



(33) 



e 



1 



4 7 2 (X F + l/2) . 



' ' " h 1 cosh(27r 7 ) - cos(2ttX f ) [ l + 4 7 2 
1 



7T 



■ arctan (coth(-7T7) tan(-7rX F )) — Int I X F 



sinh(27T7) — 7sin(27rX F ) 
1\ 1 



-a xy (Xv,i). 



(34) 
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Figure 2. The diagonal [Eq. (135 H and the off-diagonal [Eq. (136 H components of 
the conductivity tensor. The horizontal axis is the inverse magnetic field. 



As far as we know, an explicit analytic formula for a xy has never been reported thus 
far. 

In Fig. [51 we show the diagonal a xx and off-diagonal a xy components of the 
normalized conductivity tensor calculated by Eqs. (|25[) and ([3~4"| [or cquivalently, by 
Eqs. ([55]) and ([55)) below], respectively. The parameters are selected to be typical 
values in a GaAs/AlGaAs 2DES: m* = 0.067too with mo the bare electron mass, 
£f = 7.5 meV, and T = h/(2r) = 0.12 mcV. The traces basically reproduce well- 
known behavior of a 2DES in the magnetic field: the staircase with plateaus at integer 
multiples of e 2 /h for a xy and peaks at inter-plateau transition for a xx . The non- 
monotonic l/B dependence observed in a xy for B < 1 T (the depression in — a xy that 
occurs in step with the peak in a xx ) is usually not seen in the experimental traces 
for a high- mobility GaAs/AlGaAs 2DES, but can be seen in early experiments on 
Si-MOSFET [29| and is likely to be related to the short-range nature of the impurity 
potential. (See Fig.2]below for comparison with the result in the long-range potential.) 

For brevity and for the convenience in later use, we rewrite a xx and a xy in concise 
formulas, 



s (*f,7) = 



20 



1 + 4 7 2 



A F + ~ ) FmhImAV.-.) 



(35) 



a xy (X F) 7) = -IFsinh(A F , 7) - 7Fsin(A F , 7) 



4 7 2 



1 + 4 7 2 



A,, 



- )Fsinh(A F , 7 ), 



where we introduced the notations Fsin(Ap,7), Fsinh(Ap,7), 
defined as 

sin(27rAF) 



(36) 

and IFsinh(AF,7) 



Fsin(Ap, 7) 
Fsinh(A F ,7) 
IFsinh(A F ,7) = 



cosh(27T7) — cos(27rA F ) ' 

sinh(27T7) 
cosh(27T7) — cos(27rA F ) ' 

dAFsinh(A, 7) 



— arctan (coth(7T7) tan(7rA F )) + Int [ A F + — J + — . 
7r V 2 / 2 



(37) 



Although it appears, at first glance, that the stepwise behavior of a xy is reflecting 
only the first term in Eq. (|36|) , the second term is also playing its own share of roles 
by extending the width of the plateau and thus making the slope of the inter-plateau 
region much steeper than it would be were it not for the term. The steepness of the 
slope is of paramount importance in our theory that attempts to explain the behavior 
of da xy /dB. 

We note in passing that the DOS given by Eq. © can also be rewritten, following 
the same procedure as in the derivation of Eq. (|35[> . as 

D(e) = DoFsinh(X, 7 ), (38) 

where D — m* / (2irh 2 ) represents the DOS of a 2DES in the absence of the magnetic 
field, and X = e/(hio c ) — 1/2 as defined earlier. Accordingly, the cumulative number 
of states N(e) below e reads 



N(e) 



£ D(e')de' = ^IFsinh(A, 7 ). 



(39) 



We will show in |Appcndix A| that Eqs. (|35|) and (l36|) tends to the well-known 
semiclassical formulas for B — ► 0. 



4. The relation between diagonal and off-diagonal conductivities at high 
magnetic fields 



We now move on to the main topic of the present paper, the relation between a xx and 
a xy at high magnetic fields. Since both of X F and 7 are functions of B, the derivative 
of the off-diagonal component a xy with respect to B is written as 

da xy (X F ,j) _ da xy {X F ,j) dX F da xy {X F , 7) d7 



dB 



dX* 



dB 



dB 



1 

B 



1\ da xy {X Fll ) dd xy {X Fll ) 
F+ 2) dX F +7 9 7 



(40) 



Differentiation by Xp and by 7 can be analytically done on Eq. (|36[) and we obtain 



B 



da xy (X F ,j) 
dB 



X F + \ 



l-4 7 2 



1 



(l+4 7 2)2 

8 7 2 



8 7 2 

— j27T7 coth(27T7) 



2 J 1 + 4 7 2 
4 / 1 



l + 4 7 2 \ XF • 2 



I + 47 

27T7Fsinh 2 (X F ,7) 
2 



Fsinh(A F ,7) 



27r7 2 Fsin(X F ,7)Fsinh(A F ,7), (41) 



or, with the aid of Eq. ([35 



B 



da xy (X F ,j) 
dB 



1 

2^ 



l-4 7 2 



- 2tt7(1 + 87 2 )coth(27T7) 
tt (l + 87 2 )(l + 4 7 2 )^ 2 



1 + 4 7 2 



&xx{X F , 7) 



27 



X F + 1/2 



tt(1 + 4 7 2 ) 
2 sinh(2vr7) 



l + 4 7 2 



(X F + 1/2Y 



-4 



sin(27rX F )aL(X F ,7). (42) 
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We will pick out the dominant term at high magnetic fields from the r.h.s. of Eq. 
Since 7 = Y/(huj c ) tends to zero with the increase of the magnetic field, we 
expand the coefficients in terms of 7 for this purpose as 



B 



da xy (X F ^) 
dB 



2 JL 2 
3 



27 (x F 



• o( 7 ) 



4(X F 



- 1 



+ 0( 7 )}sin(2TrX F )a 2 xx (X F , 7 ). (43) 



The diagonal component a xx can be readily seen from Eq. (|35p to take peaks at 
X F = N (integer), namely when the Fermi energy lies at the center of Nth Landau 
level, with the peak height given by 



:(JV,7) = 



27 



1 



4 7 2 



X 



Ft" 



sinh(27T7) 
cosh(27T7) — 1 
8 2tt\ „ 

— + t h 

7T 3 / 



0(7 4 ) 



(44) 



and ~ away from the sharp peaks (see also Fig. [J). From Eqs. (|43[) and 
and noting that sin(27rAF) ~ at Xp ~ iV, we find that the second term in Eq. (|4"3"| 
makes the dominant contribution, leading to our final result, 

da xy (X F ,j) huj r 



dB 



irn a xx (X F ,j), 

e F 



(45) 



or A = (hj e 2 )irheij,(m* £f) _1 in Eq. Here we made use of the mobility /i = 

er/m* = eS/(2mT). Plots of da xy /dB calculated using Eq. (f4"l"T) (solid red line) and 
TTfJ,(huj c /e F )a xx with (7 xa; computed by Eq. ([55]) (dashed green line) shown in Fig. [3] 
attest to the validity of Eq. (I45|) for B > 1 T. The deviation seen at lower magnetic 
fields is attributable to higher order terms in 7 neglected in Eq. (j45|) . In Fig. we 
used the same parameter values as in Fig. [2j 

It is interesting to point out that we obtain the relation da xy /dB oc a xx instead 
of Eq. (|4"5"|) if we keep only the first term in Eq. 

da xy (X Fll ) (1) 



dB 



1 
B 



X F + ^) Fshjhi.V,-.") 



7Fsin(X F ,7) 



(46) 



which is not legitimate as discussed below Eq. (|37|) in Sec. [3l In fact, the peaks 
calculated by Eq. (|46|) exhibit much larger width and smaller (roughly half) height 
compared with those calculated by Eq. (|4"Tj) , as displayed in Fig. [3] 



5. Comparison with experimental results 

5.1. Modification for long-range potential 

In this section, we make an attempt to compare the relation between a xx and a xy 
deduced in Sec. 0] to the experimental results obtained in a GaAs/AlGaAs 2DES. 
It is well known that the main source of scattering in a GaAs/AlGaAs 2DES is the 
ionized donors. The donors are set back from the 2DES plane by a spacer layer 
with the thickness typically a few tens of nanometers. Therefore, the scattering in a 
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Figure 3. The plots of da xy /dB calculated by an unabridged equation, Eq. 14H 
[thin solid red line, plotted in both (a) and (b)] , and by an approximated equation, 
Eq. H45I I, with a xx calculated by Eq. J35P [thick dashed green line in (a)]. We also 
plot da xy /di?' 1 ' in Eq. I I46I I. obtained by keeping only the first term in Eq. (136 b . 
for comparison [dot-dashed blue line in (b)]. The traces are separately plotted in 
(a) and/or (b) for clarity. 



GaAs/AlGaAs 2DES should be described by a long-range impurity potential. Since 
a short-range potential is assumed in our theory, slight modification is necessary to 
implement the comparison. This is done by following the prescription given in Ref. 

First we observe that the off-diagonal component Eq. ([55)) can be rewritten in 
the form presented in Ref. [3], a xy — — e[dN(£-p)/dB] — uj c t(j XXi which reads in the 
normalized form, 

hdN(e F ) 1 . /A _. 
^ = —e-dB--^ axx - (47) 

The equivalence of Eq. (|47j) to Eq. (|36|) can readily be verified by performing the 
differentiation by B on N(ep) given by Eq. (|39|) : 

h dN{e F ) = IFsinh ^ p; 7 j + 7 Fsin(X F , 7) - ( X F + -] Fsinh(X F , 7). (48) 



e oB \ 2 

In a long-range potential, it is important to recall that the scattering is 
characterized by two distinct scattering times, namely, the quantum scattering time 
T q =h/(2T) that describes the impurity broadening of the Landau levels and the 
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momentum relaxation time r m = a m* / (n e e 2 ) related to the conductivity at B = 0. 
The latter time is typically 10 times larger than the former in a GaAs/AlGaAs 2DES, 



while the relaxation times are simply t„ 



t for short-range scatterers. Coleridge 



et al. [30] suggested an appropriate way of replacing t by either of r q or r m , with which 
the resultant a xx and a xy describe the conductivities under the long-range potential 
quite well. The method, in our notation, is to replace 7 only in the prefactor of Eq. 
(05]) by 7 m = 1/ (2w c r m ), leaving 7 = j q = T/ (hw c ) = l/(2uj c T q ) in Fsinh(A, 7) intact: 

27™ ( 1" 



-LR 

® XX 



(X F ,7g,7m) 



1 + 4 7 2 



X F 



Fsinh(X F ,7g) 



(49) 



The Hall conductivity is obtained by substituting 7 m and a xx ^ into the second term 
of Eq. g7D as 



~ LR 



(A F ,7g,7m) = 



h dN(e F ) 



dB 



1 



27, 



■(T LR 

XX 



= -IFsinh(A F , 7,) - 7,Fsin(X F , 7,) + 27 m CT ; 
With these substitutions, the derivative of a x ^ by B reads 

— = 7/i(A F ,7«) + V2{X F ,Jq) +%(A F ,7g,7 m ) 

with 



-LR 
xa; ■ 



(50) 



(51) 



1 



Fsinh(A F , 7,) - 7,Fsin(A F , 7,) 



V2(X F ,j q ) = i I (a f + i ) 27r7 g Fsinh(A F ,7 9 ) [Fsinh(A F , 7,) - coth(27T7 9/ , 



7 g Fsin(A F ,7 9 ) - 27T7 2 Fsin(X F , 7 9 )Finh(A F , 7^ 



and 

%(^F,7g,7m) = - 



27m - LR/' v \ 
—q-Vxx (^F,7 ? ,7m) 



x ^ 27T7gCoth(27r7g) + 1 + 



1 + 4 7n 



- 2tt 



X F + - Fsin(A F , 7) + 7 9 Fsinh(A F , 7? ) 



where the terms r/i, r\i and 773 are derived from the first, the second and the third 
term in Eq. ([50]) . respectively. Accordingly, the dominant term at high magnetic field 
changes from Eq. (|4"5)) to 



d^(A F , 74,7m) Mm^. 



dB 



\x q e F 



(52) 



[A = {h/ e 2 )TThe(ii 2 m l [i q )(m*s F )~ 1 in Eq. @], where fj, q — eT q /m* and (i m — er m jm* 
are mobilities corresponding to r q and r TO , respectively. 

In Fig.QJ we show the longitudinal and the Hall conductivities calculated by Eqs. 
(|49j) and (|50|) with parameters e F = 7.5 meV, fi q = 7.1 m 2 /(Vs) (corresponding to T 
= 0.12 meV), and \i m = 78 m 2 /(Vs). The parameters are taken from our experiment 
to be presented below. The diagonal component a xx has become much smaller than 
in Fig. [2] (note the 10 times magnification in Fig. |4j, in accordance with experiments 
in a GaAs/AlGaAs 2DES. Note that the non-monotonic behavior of a xy observed in 
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0.5 1.0 1.5 



1/5(1"') 

Figure 4. The diagonal [Eq. ]49l l1 and the off-diagonal [Eq. (1 50 M components 
of the conductivity tensor modified to account for the long-range potential. The 
horizontal axis is the inverse magnetic field. 




0.2 0.4 0.6 0.8 1.0 1.2 



1/B(T~') 



Figure 5. The plots of da^y/dB without approximation, Eq. {51} (thin solid 
red line), and da^^/dB approximated by Eq. 5 2 tl with calculated by Eq. 
J4Q 6 (thick dashed green line). 



Fig. [5] has vanished in Fig. [¥] The high accuracy of the approximation given by Eq. 
([52]) at high enough magnetic fields (B > 1 T) is demonstrated in Fig. [5] 

Although Eq. (p)2")) as well as Eq. is intended for the use in high magnetic 
fields, stringent comparison with experimental results is possible only in a rather low 
magnetic-field range (B <0.5 T) for a couple of reasons to be discussed in Sec.[6l Above 
all, we neglected the spin of the electrons altogether in the theory. In principle, the spin 
can be included in the theory by adding agfisB to En with a = ±1/2 representing the 
spin and /ib the Bohr magneton. Difficulty arises, however, because of the dependence 



of the g-factor on the magnetic field owing to the exchange interaction 31|: the g- 
factor experiences strong enhancement at the magnetic field where the Fermi energy 
lies between the Zeeman gap (exchange enhancement), which defies simple analytical 
treatment. If we limit ourselves to B <0.5 T, spin splitting can be completely neglected 
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200r 




1/B(T _1 ) 



Figure 6. The plots of da^/dB without approximation, Eq. I|51| l [thin solid 
red line, plotted in both (a) and (b)], and da]^/dB approximated by Eq. J53H 
[thick dashed green line in (a)] or by Eq. (1 52 I t [thin dotted green line in (b)] for 
a lower magnetic field range than in Fig. [5] We also plot d& xy /dB l - 1 '> in Eq. B6t 
for comparison [dot-dashed blue line in (b)]. The traces are separately plotted in 
(a) and/or (b) for clarity. 



because of the small (bare) <?-factor g = —0.44 in GaAs. In this low magnetic field 
range, the approximation in Eq. (|52p that retains only the leading term in 7 turns 
out to be insufficient, as demonstrated in Fig. [5] The approximation is improved by 
keeping the terms deriving from the first two terms in the r.h.s. of Eq. (|50p. 771 and 772, 
except for the term including Fsin(Xp,7) (the third term 773 can safely be neglected 
since j m < 7,): 



doffi(^F,7g,7m) 
dB 



(J, hhJ c r~LR, y o2 

n [<J XX (A F ,7,,7 m j 



Mr, 



1 - 



■ coth 



~LR 

® xx 



(X F ,7«.7m)- 



(53) 



fXqB ~ 

In Fig. [SI we plot da x J ^(XF,j)/dB without approximation, Eq. (|5Tj) , along with 
approximated traces, Eqs. and (|5"3"|) . calculated using ct^(Xf,7) in Eq. AS)) . 
Deviation of Eq. (|52|) from the exact result becomes evident below ^0.5 T, while Eq. 
(|53"j) reproduces the trace almost indistinguishable from that of the exact calculation 
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in the magnetic field range shown in Fig. O 

5.2. Relation between experimentally observed longitudinal and Hall conductivities 

Let us now turn to our experimental data. We prepared a GaAs/AlGaAs 2DES sample 
with fi m = 77 m 2 /(Vs) and n e = 2.1 xlO 15 m~ 2 , hence £f = 7.5 meV, shaped in a 
Hall bar geometry by photolithography. The quantum mobility fj, q = 7.1 m 2 /(Vs) 
was determined from the damping of the amplitudes ApsdH of the Shubnikov-de 
Haas (SdH) oscillation at low ma gnet ic fields, Apsda(B) / po = C expj— %/(fi q B)] 
with po the resistivity at B = [32]. It was pointed out in Ref. [32j that the 
prefactor C equals 4 in a homogeneous 2DES and the deviation from the value 
is attributable to the inhomogeneity. We have verified that the SdH amplitudes 
in our sample were described by the above equation with C = 4 reasonably well, 
confirming that inhomogeneity is minimal in our sample. (Note, however, that small 
inhomogeneity is inevitably present in a 2DES grown by molecular beam epitaxy, as 
will be discussed below.) Measurements were done in a dilution refrigerator equipped 
with a superconducting magnet at the base temperature (~15 mK), a temperature 
low enough for the approximation k^T <C £f to be valid. The standard low-frequency 
(13 Hz) ac lock-in technique was employed for the resistivity measurement with a low 
excitation current (10 nA for B <^ 1 T and 0.5 nA for higher magnetic fields) to 
prevent the electron heating. For the magnetic field sweep, we adopted very slow 
sweep rates (0.01 T/min for B <~ 1 T and 0.1 T/min for higher magnetic fields), 
which, combined with a high data acquisition rate (~ 4 data points/s), allow us to 
acquire data points dense enough to perform the numerical differentiation with respect 
to B reliably. The slow sweep rates are also favorable in avoiding the hysteresis in 
the superconducting magnet that obscures the exact value of the magnetic field felt 
by the sample. The longitudinal and the Hall resistances measured in our Hall bar 
sample are translated to resistivities p xx and p xy by using the geometrical factors of 
the Hall bar. Then we obtained o~ xx and o~ xy by numerically inverting the tensor, 
o-xx = Pxx/{p 2 xx + Ply) and o- xy = Pxy/{p 2 xx + Pxy)- As mentioned earlier, spin 
splitting can completely be neglected for B < 0.5 T. Due to the spin degeneracy, the 
conductivities experimentally measured in this magnetic field range are simply twice 
as large as those without the spins; considering the spin degeneracy, the normalized 
conductivities are defined here as a a p — a a p/ (2e 2 /h). 

In Fig. [7J we show da xy /dB attained by the numerical differentiation of 
experimentally obtained o~ xy , and da xy /dB approximated by Eqs. (|52| and (|53| using 
experimentally acquired o~ xx . It can be seen by comparing Figs.[H]and[7]that our theory 
reproduces the experimentally obtained traces remarkably well. Note that the same 
vertical scale is used for the two figures. Both figures reveal that the approximation by 
Eq. ((52]) progressively worsens with decreasing magnetic field, while Eq. ([53]) remains 
a good approximation over the magnetic field range shown in the figure. We want 
to emphasize that the good quantitative agreement, demonstrated in Fig. [TJ between 
da xy /dB directly deduced from o~ xy and that approximated by Eq. (|53[) using o~ xx is 
achieved without any fitting parameter. In Fig. [3 we also plot the r.h.s. of Eq. ([3]), 
a more conventional empirical relation. For the coefficient /3, we adopted the relation 
(3 = 2r m /r 9 = 2p, m /p q proposed by Coleridge et al. [5j. We can see that Eq. ([53]) 
describes the relation between o~ xx and a xy much better than Eq. ([3|) . It is clear from 
the figure that even if we use (3 as a fitting parameter, agreement by Eq. ^ cannot 
be improved very much. 
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Figure 7. Experimental traces to be compared with Fig. [6] dcr xy /dB deduced 
by numerical differentiation of experimentally obtained cr xy [thin solid red line, 
plotted in both (a) and (b)], da xy /dB approximated by Eq. J53D [thick dashed 
green line in (a)] or by Eq. (1 52 I t [thin dotted green line in (b)] calculated using 
experimentally obtained a xx . The r.h.s. of Eq. Q with the experimentally 
obtained cr xx and (3 = 2Tm/r q is also plotted by dot-dashed blue line in (b). 
The traces are separately plotted in (a) and/or (b) for clarity. 

For higher magnetic fields, spin splitting manifests itself as the splitting of the 
peaks in a xx and da xy /dB. The peaks take place at the conditions £f = EN+g*<J[iBB 
[N = 0, 1, 2,..., a = ±1/2, and g* represents the g-factor including (independent) 
exchange enhancement], instead of ep = En in the spin-degenerate case, and therefore 
Eqs. (|49|) and (f50|) no longer describe the positions of peaks or steps between adjacent 
plateaus correctly. Nevertheless, concurrent occurrence of peaks in a xx and in da xy /dB 
still allows us an attempt to see the applicability of Eq. ([52"]) . as shown in Fig. [5] 
Here d a p — <r a f]/(e 2 /h) again since spin degeneracy is now lifted. We see that Eq. 
(p)2"j) reproduces roughly the right order of magnitude for the height of the peaks in 
da xy /dB, although the increase in the peak height with increasing magnetic held for 1 
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Figure 8. Experimental traces for higher magnetic-field range: da xy /dB deduced 
by numerical differentiation of experimentally obtained a xy [thin solid red line, 
plotted in both (a) and (b)], da xy /dB approximated by Eq. H52H calculated using 
experimentally obtained a xx [thick dashed green line in (a)], and the r.h.s. of Eq. 
(O with the experimentally obtained cr xx and /3 ~ 400 [dot-dashed blue line in 
(b)]. The traces are separately plotted in (a) and/or (b) for clarity. 



T<_B<2.5Tisat obvious variance with the behavior of da xy /dB. The discrepancy 
is mainly ascribable to the deviation of experimental peak heights in a xx from the 
ocl/B dependence inferred from Eq. (I49|) . By contrast, we find that our experimental 
result is well described by Eq. ([3|) in accordance with previous studies [2,[3|,U,[5|,[a,LZl,l8|], 
albeit with the value of the parameter (3 ~ 400 roughly 20 times larger than 2r TO /r g . 

6. Discussion 

The relation between a xx and a xy is already implicit in Eqs. (|18[) and (|19[) , since 
both of the components derive from the same set of Green's function Gn(£) (N = 
0, 1, 2,...), or from the same DOS; note that once the DOS is given, both imaginary 
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and real parts of Gn(s) are known by Eq. pip and by the Kramers-Kronig relation, 
respectively. We have shown in Sec. Q] that the relation can be explicitly written 
down as Eq. (|42p . or approximately as Eq. (145p . if we assume a simple form given by 
Eq. (|12p for the Green's function corresponding to the Lorentzian broadening of the 
Landau levels Eq. flSJ. It might be argued that the expression Eq. (fl~2"|) is too crude 
to represent a 2DES under magnetic field. We expect, however, that improvement in 
Gn(s) does not alter the relation Eq. (|42|) to a large extent [if we keep ourselves within 
the framework of the approximate relation between the conductivity tensor and the 
Green's function represented by Eqs. (|18p and (|19p ], so long as the resultant DOS does 
not significantly deviate from the Lorentzian line shape. An important point we would 
like to stress is that the relation Eq. (|42|) is inherent in the expressions of a xx and a xy 
and requires no external source, e.g., the inhomogeneity in the electron density. 

In Eqs. (fP8|) and (fl"9|) . we have neglected a number of effects known to take place 
in a 2DES subjected to a magnetic field. These include the localization, the formation 
of the edge states (stripes of compressible states parallel to the edge of the sample 
interleaved with incompressible regions), and the electron-electron interaction. We 
have also neglected spins altogether as mentioned in Sec. [5l Due to the localization in 
the tails of Landau level peaks, the width of the peaks in a xx will become narrower than 
what is shown in Figs. [5] and U for the high magnetic field region where overlap between 
adjacent Landau level peaks can be neglected. The electron-electron interaction will 
engrave additional minima on the peaks of a xx between adjacent integral quantum 



Hall states for N < 2 Landau levels via the fractional quantum Hall effect [331 ] . 
and also affect the height and shape of the peaks for higher Landau levels through 
forming the (probably incomplete) charge density wave states IH 35| . For long-range 
impurity potential, the peaks will be altered also by the network of the compressible 
and incompressible stripes formed around valleys or hills of the impurity potential 
[36j |. Strictly speaking, therefore, our theory applies only to the low-magnetic field 
region where these effects are negligibly small. This is exactly the region we have 
employed in the comparison with the experimental result in Fig. [7J The excellent 
agreement between the theory and experiment attests to the correctness of our theory 
were it not for the additional effects neglected in the theory. The slight difference 
in the line shape between theoretical (Fig. [6]) and experimental (Fig. [7]) traces, with 
the theoretical trace showing asymmetry between sharp maxima and rather rounded 
minima, is attributable to the use of constant £f in the theory; in the experiment, £f 
is expected to oscillate with magnetic field to keep the electron density n e constant, 
resulting in more symmetric peaks and dips [37| . 

In the higher magnetic-field regime, we envisage better agreement between 
theoretical and experimental results by modifying our theory to include the effects 
neglected in the present paper listed above, which is the subject of our future study. 
In the high-magnetic-field regime, however, we are unable to rule out the possibility 
that the inhomogeneity in n e is the dominant source of the experimentally observed 
relation Eq. (fTJ) [or Eq. (J3]) as shown in Fig. [8], as suggested by previous studies 
@> HI 11 1; the effect of the inhomogeneity is expected to gain more significance at 
higher magnetic fields, since the difference in the Hall resistivity Ap xy between two 
points differing in the electron density by An e , Ap xy ~ An e B / (n^e), increases with 
B. Note that, in realistic samples, both microscopic inhomogeneity owing to the 
random distribution of the dopants and macroscopic inhomogeneity resulting from 
the technical difficulties in the molecular beam epitaxy are virtually impossible to be 
completely eliminated. 
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7. Conclusions 



We have calculated the diagonal (<r xx ) and off-diagonal (a X y) components of the 
conductivity tensor in the quantum Hall system by the linear response theory, 
neglecting the correction from the current vertex part. A Lorentzian line shape with 
the width Y independent of the magnetic field was assumed for the broadening of 
the Landau levels by the short-range impurity potential. The corresponding simple 
approximation for the Green's function Eq. (|12p allowed us to obtain analytic formulas 
for both a xx and <r xy , given by Eqs. (|3"5")l and (|3"o]l respectively, for k^T <C £f- The 
formulas asymptotically approach the semiclassical formulas at low-magnetic fields. 
Inspection of the formulas reveals that da xy /dB is proportional to Ba xx [Eq. (|45p] at 
high magnetic fields where Y <C fiu> c - This comprises a possible alternative route to 
explain, without resorting to the inhomogeneity in the electron density, the well-known 
empirical relation between a xx and a xy . 

To account for the long-range nature of the impurity potential in a GaAs/AlGaAs 
2DES, slight modification was made by introducing two types of scattering times, the 
quantum scattering time r q and the momentum relaxation time r m , yielding Eqs. 
(|49p and ([50)1 for a xx and a xy , respectively. The resultant relation between the two 
components, Eq. (|53j) . is found to be in quantitative agreement with the experimental 
result obtained in the GaAs/AlGaAs 2DES at the magnetic field range where the 
spin splitting, the localization, the formation of the edge states, the electron-electron 
interaction, etc., can be neglected. 
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Appendix A. Conductivity tensor in the weak magnetic field limit 



In the derivation of Eqs. (|35|) and (J36J), we have made no assumption on the 
strength of the magnetic field. Therefore the equations should, in the low-field limit, 
asymptotically coincide with the well-known semiclassical expressions 

°™ = 1 + (<lr)2 (A - 1} 



sc _ a uj c T n e e 1 sc 



with do = n e e rjm* — spe t / '(2 / !th , ) or, in the normalized forms, 



~sc 2 7 ( v , 1 



^ = irtf i Xf + - ' (A - 3) 
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3% = (*f + ~) = - (x F + ~) + 2 7 S*° (A.4) 

Since Fsinh(AF, 7) — ► 1 with 7 — * 00, it is ready to see a xx — > with £> — * in Eq. 
([35]) . From Eqs. ® and (36]), we find 

c^,(X F , 7 ) = -IFsinh(X F ,7) - 7Fsin(A F , 7 ) + 2-fa xx (X F ,j). (A.5) 

Noting that IFsinh(A F , 7) — > (X F + 1/2) and 7Fsin(A F , 7) — > with 7 — > 00, we can 
also perceive a xy — > a X y with B — > 0. 
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